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H . Abstract 

In this work we investigate the background dynamics when dark energy is coupled to dark matter 
with a suitable interaction in the universe described by brane cosmology. Here DGP and the RSII brane 
a^: models have been considered separately. Dark energy in the form of modified Chaplygin gas is considered. 

A suitable interaction between dark energy and dark matter is considered in order to at least alleviate (if 
O ■ not solve) the cosmic coincidence problem. The dynamical system of equations is solved numerically and 

a stable scaling solution is obtained. A significant attempt towards the solution of the cosmic coincidence 
5h . problem is taken. The statefinder parameters are also calculated to classify the dark energy models. 

Graphs and phase diagrams are drawn to study the variations of these parameters. It is also seen that the 
background dynamics of modified Chaplygin gas is completely consistent with the notion of an accelerated 

■ expansion in the late universe. Finally, it has been shown that the universe in both the models follows 
^ , the power law form of expansion around the critical point, which is consistent with the known results. 

oo 

^ ! 1 Introduction 

The idea of the static universe once conceived by Albert Einstein faced a big threat when at the turn of the last century 

■ observations of la Supernova confirmed that our universe is suffering from an accelerated expansion (Perlmutter, S. 
et al., 1999; Spergel, D. N. et al.2003). Now belief and observations are always to be supported by some model. In 
the quest of finding a suitable model for universe, Cosmologists started to investigate the root cause that is triggering 
this expansion. Fundamentally, we were to modify Einstein's equation either by modifying the left hand side, i.e., 

■ modifying the idea of Einstein gravity or to modify the right hand side which immediately speculate the nature of 
5-H , the matter inside the universe. If our Universe is filled up by some invisible fluid causing a negative pressure then it 
^ ■ violates the strong energy condition i.e. p + 3p < 0. Because of its invisible nature this energy component was aptly 

termed as dark energy (DE) (Riess, A. G. et al., 2004). 

With the introduction of DE, the search began for different candidates that can effectively play the role of DE. 
DE represented by a scalar field □ is often called quintessence. Not only scalar field but also there are other Dark 
fluid models like Chaplygin gas which plays the role of DE very efficiently. The earliest form of this was known as 
pure Chaplygin gas (Kamenshchik, A. et al.,2001; Gorini, V. et al. 2004). Extensive research saw pure Chaplygin 
gas first modify into generalized Chaplygin gas (Gorini, V. et al., 2006; Alam, U. et al. 2003; Bento, M. C. et al. 
2002;Carturan, D. and Finelh, F. 2003; Barreiro, T. and Sen, A.A. 2004) and later into modified Chaplygin gas (MCG) 
obeying the equation of state (Benaoum, H. B., 2002; Debnath, U., 2004) 

P = Ap-^, (1) 
P 

where p and p are respectively the pressure and energy density and < a < \, A and B are positive constants. 
Currently, we live in a special epoch where the densities of DE and DM are comparable. Although they have evolved 
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independently from different mass scales. This is known as the famous cosmic coincidence problem, till date several 
attempts have been made to find a solution to this problem (del Campo, S. et al 2009; Leon, G. and Saridakis, E.N. 
2010; Jimenez, J.B. and Maroto, A.L.2008; Berger, M.S. and Shojae, H. 2006; Zhang, X. 2005; Griest, K. 2002; Jamil, 
M. and Rehman, F. 2009; Jamil, M. et al. 2010a; Jamil, M. and Saridakis, E. N. 2010b; Jamil, M. and Farooq, M. 
U., 2010c; Jamil, M. et al. 2010d). A suitable interaction between DE and DM is required if we wish to find an 
effective solution to this problem. It is obvious that there has been a transition from a matter dominated universe 
to dark energy dominated universe, by exchange of energy at an appropriate rate. Now the expansion history of the 
universe as determined by the supernovae and CMB data (Jamil, M.et al. 2010a; Jamil, M. and Saridakis, E. N. 
2010b) bounds us to fix the decay rate such that it is proportional to the present day Hubble parameter. Keeping the 
fact in mind cosmologists all over the world have studied and proposed a variety of interacting DE models (Setare, M. 
R. 2006,2007c; Hu, B. and Ling, Y. 2006; Wu, P. and Yu, H. 2007; Jamil, M. 2010e, 2011; Dalai, N. 2001). 

Here we shall use MCG as the DE and the equation of state for MCG is given by ^ with p = pmcg Recent 
developments regarding MCG can be found in (Lu, J. et al. 2008; Dao-Jun, L., Xin-Zhou, L. 2005; Jing, H. et. al. 
2008; Debnath, U. 2007). This special form also appears to be consistent with the WMAP 5-year data and henceforth 
support the unified model with DE and DM based on generalized CG (Makler, M. et al. 2003; Setare, M.R. 2007a, 
2007b, 2009; Barreiro, T. et al. 2008). Recent supernovae data also favours the two-fiuid cosmological model with CG 
and matter (Panotopoulos, G. 2008). 

As we have stated earlier, modifying the right hand side of Einstein's equation was not the only way to explain 
the increase in the rate of the expansion, rather to modify the gravity part of the left hand side is also to demonstrate 
the present day universe. In this branch Brane-gravity was introduced and brane cosmology was developed. A review 
on brane-gravity and its various applications with special attention to cosmology is available in (Rubakov, V. A. 2001; 
Maartens, R. 2004; Brax. P. et al. 2004; Csa ki, C. 2004). In this work we consider the two most popular brane 
models, namely DGP and RS II branes. Our main aim of this work is to examine the nature of the different physical 
parameters for the universe around the stable critical points in two brane world models in presence of MCG. Impact 
of any future singularity caused by the DE in brane world models will be studied. 

This paper is organized as follows: Section 2 comprises of the dynamical system analysis in DGP brane model. 
Section 3 deals with the dynamical system analysis in RS II brane model, followed by a general discussion in section 
4. 

2 Dynamical system analysis in DGP brane world model 

A simple and effective model of brane-gravity is the Dvali-Gabadadze-Porrati (DGP) braneworld model (Dvali, G. 
R. et al. 2000; Deffayet, D. 2001; Deffayet, D. et al. 2002) which models our 4-dimensional world as a FRW brane 
embedded in a 5-dimensional Minkowski bulk. It explains the origin of DE as the gravity on the brane leaking to the 
bulk at large scale. On the 4-dimensional brane the action of gravity is proportional to whereas in the bulk it is 
proportional to the corresponding quantity in 5-dimensions. The model is then characterized by a cross over length 

scale Tc = -^W ^^^^ that gravity is 4-dimensional theory at scales a «rc where matter behaves as pressureless dust, 
but gravity leaks out into the bulk at scales a» Tc and matter approaches the behaviour of a cosmological constant. 
Moreover it has been shown that the standard Friedmann cosmology can be firmly embedded in DGP brane. 

It may be noted that in literature, standard DGP model has been generalized to (i) LDGP model by adding a 
cosmological constant (Lue,A. and Starkman,G. D. 2004), (ii) QDGP model by adding a quintessence perfect fluid 
(Chimento,L. P. et al. 2006), (iii) CDGP model by Chaplygin gas (Bouhmadi-Lopez, M. and Lazkoz, R. 2007) and 
(iv) SDGP by a scalar field (Zhang, H. and Zhu, Z. H. 2007). In (Wu, X. et al. 2008) the DGP model has been 
analysed by adding Holographic DE (HDE). 

2.1 Basic equations in DGP brane model 

While flat, homogeneous and isotropic brane is being considered, the Friedmann equation in DGP brane model(Dvali, 
G. R. et al, 2000; Deffayet, D., 2001; Deffayet, D. et al, 2002) is modifled to the equation 



where H — the Hubble parameter, p is the total cosmic fluid energy density and Tc — i^jfi is the cross-over scale 
which determines the transition from 4D to 5D behaviour and e = ±1 (choosing — SttG — 1). For e — +1, we 
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have standard DGP(+) model which is self accelerating model without any form of DE, and effective w is always 
non-phantom. However for e — —1, we have DGP(— ) model which does not self accelerate but requires DE on the 
brane. It experiences 5D gravitational modifications to its dynamics which effectively screen DE. Brane world scenario 
is actually a modified gravity theory. If we write the Einstein equation for brane world in terms of Einstein gravity 
then the extra term can be treated as the effective DE But that is not the physical DE. Moreover this DE is applicable 
only in Einstein gravity. But here we will consider the physical DE in brane world. So we have introduced the MCG 
type fluid in brane. There are few many candidates of DE models, such as scalar field (Nojiri, S.and Oditsov, S.D. 
2004), phantom (Nojiri, S. et al 2005), k-essence (Bamba, K.et al 2011), tachyonic field ((Nojiri, S.and Oditsov, S.D. 
2003), etc. MCG being a modified form of DE with a unique equation of state for DE and DM. Also MCG is the 
observational viable model (Thakur, P. et al 2009, Dev, A. et al 2003). It interpolates between radiation in early stage 
and ACDM in late stage. 

As in the present problem the interaction between DE and pressureless dark matter has been taken into account 
for interacting DE and DM the energy balance equation will be 

Pmcg +3H {1+ Wmcg) Prncg -Q, for MCG and (3) 

Prn + SHp„i = Q, for the DM interacting with MCG. (4) 

where Q — SbHp is the interaction term, b is the coupling parameter (or transfer strength) and p — pmcg + Pm is 
the total cosmic energy density which satisfies the energy conservation equation p + 3H {p + p) =0 (Guo, Z.-K. and 
Zhang, Y.-Z. 2005; delCampo, S. et. al. 2008). 

As we have lack of information about the fact how do DE and DM interact so we are not able to esteemate the 
interaction term from the first principles. However, the negativity of Q immediately implies the possibility of having 
negative DE in the early universe which is overruled by to the necessity of the second law of thermodynamics to be 
held(Alcaniz, J. S. and Lima, J. A. S. 2005). Hence Q must be positive and small. From the observational data 
of 182 Gold type la supernova samples, CMB data from the three year WMAP survey and the baryonic acoustic 
oscillations from the Sloan Digital Sky Survey, it is estimated that the coupling parameter between DM and DE must 
be a small positive value (of the order of unity), which satisfies the requirement for solving the cosmic coincidence 
problem and the second law of thermodynamics (Feng, C. et. al. 2008). Due to the underlying interaction, the 
beginning of the accelerated expansion is shifted to higher redshifts. Consequently using the Friedmann equation ([2]) 
and the conservation equation, we obtain the modified Raychaudhuri equation 

2H^-)h = -H{p + p), (5) 
2.2 Dynamical system analysis 

In this subsection we plan to analyse the dynamical system. For that firstly the conversion of the physical parameters 
into some dimensionless form, given by 

, Pmcg Pm 

a; = lna, u = ^, v = tt d 

where the present value of the scale factor oq = 1 is assumed. Inclusion of ([T]) into ([5]) gives us the parameter gradients 
as 

du 

dx ' V ' J ' e{u + v) - + (e2 - l)(u + v) 

and 



3^ {e - VI + (e^ - 1)(" + v)} {u{l + wifcT'') + v] 
3b {u + v)-3u(l + w^,ZT^ ) + ^ ^ ^ = ^ (7) 



dy SV \e-^l + (£2 - l){u + V)} \uil + W^mc^^^) + v} 

^^3b{u + v)-3v + ^ ^ ^ = ^. (8) 

dx eiu + v)- + -l)iu + v) 

(DGP) 

Where, Wmcg is the EoS parameter for MCG determined as 

(a+l) 



(DGP) ^ ^A-B 

meg -i i- 

Pmcg 



2'^rl{\- u - v) 



2,u (e - Vl + (e2 -!)(« + «))' 



(9) 
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2.2.1 Critical point 

For mathematical simplicity, we consider a = 1 only. The critical points of the above system are obtained by putting 
^ = = ^ . Now this system of equation does not yield an explicit solution. So we have to put some special values 
to the constants which will yield a non zero positive solution of the system of equations. Considering ^4 = and e = 1 
we obtain the following critical value to the system. 



Uc = 1 — b 



3D V3_ /3(1 - 36 + 362 - b^) 4 (1 + 4fe - 6b^ + 4&3 - 64) 



Br'^ 



D 



1 /9-186 + %2 3 \ , , VSTV 1 + 35 - 362 + 63 4(1-26 + 62) 6(1 + 6) 
+ -£> +3 + 6+— + — ^ ^ - 



Br-^ 



(10) 
(11) 



where 



D = (1 - 36 + 362 - 63) 



and 



,12/ , ^ 3Vl - 36 + 362 _ 53 \ 
N={ — \uc + b-l + — — ) 

The critical point correspond to the era dominated by DM and MCG type DE. For the critical point (uc,Wc), the 
equation of state parameter ^ of the interacting DE takes the form 



R92(q+1) 2(q + 1) /, \2{a+l) 
yj(DGP) ^ ^ U - Me - Vc) 

3"+iu?+i(e - ^l + (e2-l)(ue + t;e))2(a+l) 



(12) 



2.2.2 Stability around critical point 

Now we check the stability of the dynamical system (eqs. ([7]) and ([8|)) about the critical point. In order to do this, 
we linearize the governing equations about the critical point i.e.. 



u = Uc + 5u and v = Vc + Sv, 



Now if we assume f = ^ and g — then we may obtain 



and 



dv \ 



dx I 



= [du9\^5u+[d^g\Jv 



where 



(13) 

(14) 
(15) 



duf 



+ 



-36- 3Bu 



1 



4^4 



-326l4 (e2 - 1) 

9u'^K (e - Kf ^ 9^2 (e - xf ^ 9^3 (e - Kf ^ 
16B9*rt ( 32B9^ {e^ - l) 



-3 1 + 



e{u + v)- 3uK (e - K) 



1 + A- 



9u2 (e - K f \ 9u^K (e - K f 



WBO'^r, 



4^4 



MBO^rt 2,2Be^rt 



9u2 (e - if)* 9u3 (e - Kf 



3it ( e 



16-B6l*r* 



3u {e^ -l){y + u{l + A 



{{u + v)e-Ky 



2Kiiu + v)e-K) 



3{e-K)i^v + u(f + A 

{u + v) e — K 



(16) 
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-36 — 3u 



( 



(e^ - 1) 7-4 64503^4 



+ 



37i (e - A ) (^i + 7/ 9„2jf(,_K)S' + 9uHe-Kf 



{u + v) e — K 



3u 



(,_4_i)(,_i^)(, + ,(l + A- Jg^)) 37.(62-1) (. + 7.(1 + 



((u + v)e-Ky 



2K{{u + v)e-K) 



(17) 



du9 



36- 



e (7i + 7;) - 3uif (e - 



1 + 



lese'M / 32B6l4 (e2 _ 1) ^4 647*561 V + 3256*^4 ' 



3^ (^-4g^) (^--^) 



7; + 7i ( 1 + A - 



9«2(c_i<-)- 



■)) 37; (e2 - 1) (7; + w (1 + 



A 



9«2(e__R-)" 



((7t + 7;)e-iC)' 



2i4: {{u + v)e-K) 



(18) 



37> (e - K) 



-3 + 36 



1 + 77 



32Be-^(e"-l)r^ _64Befr^ 



(it + Ti) e — if 



(77 + w) e — _ftr 



37; (e - JC) 
((7t + 7;)e- X)^ 



(e2-l)W / 16B04 4 

^ ^' 7; + 77 l + yl " 



(e2 - 1) (^7; + 7i (^1 + 



9n2(e-K)* 



2a: 



97/2 (, _ xy 



2K{{u + v)t- K) 



(19) 

where K = yjl + {u + v) {e^ - 1), and 61 = (1 — 77 — 7;). The subscript c refers to quantities evaluated at 
the critical point of the dynamical system. 



2.2.3 Nature of cosmological parameters 

We calculate the deceleration parameter q = — 1 — {H/H'^), in this model as 



o I 127-2 ^ _^ 

q(DGP) ^ _ij^±) 



V4r2 + 3a UAt^ + 3a + ^ 



(20) 



where a = - . The above deceleration parameter can be written in terms of dimensionless density parameter ^mcg 



JDGP) _ 



^^3 f 12r2 (1 + 7.(^.^)17 



V4r2+3a(^V4'-?+3a+^ 



(21) 



Now using © we obtain Vlmcg = = So from equation (PT|) it is obvious that, 



12r2 1 



JDGP) 



u-\-v 



V4r2+3a(V4r2+3a+^) 



(22) 



The solution corresponding to the stable critical point, such that {u,v) — > {uc,Vc)- Hence using (|22l) we get 

„(-DGP), 



12r^ 1 



^(DGP) ^ ^ ^ ,u'/iere X(dgp) - r , . , s ■ 

2 |y4^T^(y4;:T^+_^)| 



(23) 



5 




Fig 1 : The dimensionless density parameters are plotted against e-folding time. The initial condition is 

v{0) = 0.9, u(0) = 0.2. Other parameters are fixed at a = l,b = 0, A = 1, B = 0.5 and Vc = 10. 

Fig 2 : The dimensionless density parameters are plotted against e-folding time. The initial condition is 

v(0) = 0.6, u{0) = 0.6. The other parameters are fixed at a = l,b = 0, A = 1, B = 0.5 and Vc = 10. 

Fig 3 : The phase diagram of the parameters depicting an attractor solution. The initial conditions 

chosen are v{0) = 0.5, u(0) = 0.6 (green); v{0) = 0.6,^(0) = 0.6 (blue); v{0) = 0.7,^(0) = 0.6 (red); 

v(0) = 0.8, u{0) = 0.6 (brown). Other parameters are fixed at a = 1, 6 = 0, ^ = 0.3, B = 1 and rc = 10000. 




Fig 4 : The deceleration parameter is plotted against the state parameter. Other parameters are fixed at 
B = l,b = 0,rc = 10000, e = -1. 

Fig 5 : The statefinder parameter r is plotted against the state parameter. The other parameters are fixed 
at 5 = 1, 6 = 0, rc = 10000, e = -1. 

Fig 6 : The statefinder parameter s is plotted against the state parameter. The other parameters are fixed 
at S = 1, 6 = 0, rc = 10000, e = -1. 

Fig 7 : The ratio of density parameters is shown against e-folding time. The initial conditions chosen are 
v(0)=0.6, u(0)=0.6. The other parameters are fixed at a = 1, A = l,b = 0, B = 0.5, rc = 10. 
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Moreover the Hubble parameter can be obtained as 



H= .^^ . (24) 



where we have ignored the integration constant. Integration of (j24p yields 



a{t) = aof'^'i^op^ ^ (25) 

which gives a power law form of the expansion. In order to have accelerating universe 3X(^^cp} ^ ^ ^ 
X{DGP} < §■ Using this range of Xj^jcp) in the equation q[^^^^ = — 1 + ^Xi^jjQpy We get the range of qc^^^^ as 

( DGP) 

— l<qc < 0. Therefore the deceleration parameter is negative and hence the result is consistent with the fact 
that the universe is undergoing an accelerated expansion. 

Sahni et al (2003) introduced a pair of cosmological diagnostic pair {r, s} which is known as as statefinder pa- 
rameters. The two parameters are dimensionless and are geometrical since they are derived from the cosmic scale 
factor alone. Also this pair generalizes the well-known geometrical parameters like the Hubble parameter and the 
deceleration parameter. The statefinder parameters are given by 

In the DGP brane model, we have the following expressions of r and s as 



riDGP) =[l ^ j (1 - 3X(^GP)) . (27) 

and 



S{DGP} — X(DGP)- (28) 



3 Dynamical system analysis for RS II brane world 

Randall and Sundrum (1999a, 1999b) proposed a bulk-brane model to explain the higher dimensional theory, popularly 
known as RS II brane model. According to this model we live in a four dimensional world (called 3-brane, a domain 
wall) which is embedded in a 5D space time (bulk). All matter fields are confined in the brane whereas gravity can 
only propagate in the bulk. The consistency of this brane model with the expanding universe has given popularity to 
this model of late in the field of cosmology. 



3.1 Basic equations in RS II brane model 

In RS II model the effective equations of motion on the 3-brane embedded in 5D bulk having Z2-symmetry are given 
by [Maartens, R. 2000, 2004; Randall and Sundrum, 1999b; Shiromizu et al, 2000; Maeda et al, 2000; Sasaki et al, 
2000] 

^^^C?M^' = -^iQf^i' + «:4^M"> + ^^(^•^ - Et,^ (29) 

where 

^l = \ \4 , (30) 

A4 = i«i(^A5 + i«iA2) (31) 
and 111 1 

Il^,u = -^ T^aT^" + TT^u + ^ q^,uTo.0T'^^ " ^ qt^^T^ (32) 

and E^^ is the electric part of the 5D Weyl tensor. Here K5, A5, t^^, and A4 are respectively the 5D gravitational 
coupling constant, 5D cosmological constant, the brane tension (vacuum energy), brane energy-momentum tensor and 
effective 4D cosmological constant. The explicit form of the above modified Einstein equations in flat universe are 
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2A 



and 



The dark radiation U obeys 



2 ^ 



w + mu = 



(33) 
(34) 

(35) 



where p = Pmcg + pm and p = Pmcg + Pm are the total energy density and pressure respectively. The continuity 
equations for dark energy and dark matter are given in ([3]) and (|4]). Now we shall study the dynamical system 
assuming A4 = = (in absence of cosmological constant and dark radiation). 



3.2 Dynamical system analysis 

Now to analyze the dynamical system we define three variables x, u and v as given in equation ([5]). Now making use 
of ©, ©, dlD, dSni) and ^ we get, 



§ = -3& (u + «) - 3^ (1 + +3.^^ {u (1 + u;(fj")) + .} - -i^ {l - ^^(u + «)} {u (2 + u;™) + 2.} 

(36) 



and 

dv 
dx 

where 



36 (u + w) — 3w + 3ii 



nl [u (1 + z«(„^f ^)) + - -J- {1 - .|(« + .)} (2 + «;(f//^)) + 2.} 



(RSII) ^ ^ _ 



BkI {u + ■(;) 



2uA{l - K| (u + v)) 



a+l 



(37) 



3.2.1 Critical point 

Similar to DGP brane model, the critical points of the above system are obtained by putting = = Now this 
system of equation does not yield an explicit solution. So we have to put some special values to the constants which 
will yield non-zero positive solution of the system of equations. Considering v4 = 0, i? = 1, 6 = 0.5, A = 1, a = 1 and 
K = 1 we obtain the following critical points to the system. 



uic = 0.491891, uic = 0.190983 

M2c = 1.19631, U2c = 0.836396 
U3c = 2.22826, W3c = 10.9091 



(38) 

(39) 
(40) 



3.2.2 Stability around critical point 

Now if we write / = 4^ and g = 4^, then we can obtain the following expressions 



duf 



5u 



dvf 



Sv 



and 



where 



<5 ( = [dug]cSu+[dyg]Jv 



(41) 
(42) 



duf^ 



4u2 {u + vy {{u + u)fc2 - 1}^ A^ 



[B{u + w)4fc4 {~6u^ - 2uv + v'^ + {u + v) (lOu^ - v'^) P - 4u2(u + vfk^] 
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-4:u^{{u + v)k^ - 1}^ |-3u2 - 2Au^ - bu^ - Guv - iAuv - 2buv - 3v^ ~ Av^ - hv^ + [u + vf ((6 + 4A) u + 3u) fc^j 

(43) 



B{u + vfk^ |-7m -Av + {u + v) (11m + 2v) - 4u (u + vf fe'^j 



dug = 3 



4w (u + vy {{u + v)k'^ - lyx^ 

+Au{{u + v)k^ - 1}^ i^Au^ + {u + vf {-b + 3uk^) I (44) 
3v {u + vf k'^ |-6w + 3v + {12u - 7v) {u + v)P + A{v- u) (u + vf k'^^ 



{ 3Av ^ B 
b + v\{7 + AA)k^-^^^ , , 3 



[u + vY] Au^ {{u + v)k^ -ly \^ 

(45) 



a„5 = 5 ^ ^ — (u + vf k'^ (-3 (u + 4w) + (u + w) (7ii + 2&v) fc^ - 4 (u + w)^ (u + 3w) fc'^j 

4M(u + w)'{(w + t;)fc2-l}3A2 L ^ ^ I ^ W 

+4w {(u + w) fc^ - l}^ |-7w^ - iAu^ + foti^ - 14™ + 2buv - 7v^ + bv"^ + {u + vf {(7 + AA)u + lAv} fc^j A^ (46) 



3.2.3 Nature of cosmological parameters 

In this RSII model, the deceleration parameter q can be obtained as 



3 {^Kr^^-2)-(l + ^l^r^^)} 



{B.SII) ^ _ 1 _ r V : 1 ^ V I '> (47) 

' 2 (1 + ^) 

which can be written in terms of dimensionless density parameter ^mcg = ^'p"" as in the following 

Now since ^rncg = = and assuming ^ = (-(rsu) we get, 



3{(l-^)-^'^r'^ + (l + ^(«^//))} 
-1 H ^ 



y 2 (i + f(££iii) ^ > 

Considering only the first stable critical point, such that (u, w) (uic, wic), using (I49p we get, 

If t(Rsii) = — 1^ (Rsj'j)\^ ^""^ 1 ^(RSii) — 0, we have q = —1, which confirms the accelerated expansion of the 
universe. When e^Rsii) = ^2 we have q — —oo. Therefore we have super accelerated expansion of the universe. 

In this scenario, the Hubble parameter can be obtained as, 
where the integration constant has been ignored. Integration of ([5T|) yields 



a{t) = aoi'^'t^s//) (52) 
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Fig. 8 Fig. 9 Fig. 10 

Fig 8 : The dimensionless density parameters are plotted against e-folding time. The initial condition is 
v{0) = 0.4, u{0) = 0.1. Other parameters are fixed at a = 1, 6 = 0.5, A = 0, B = 1, X = 1, and k = 1. 
Fig 9 : The dimensionless density parameters are plotted against e-folding time. The initial condition is 
v{0) = 0.3, u{0) = 0.3. The other parameters are fixed at a = 1, 6 = 0.5, A = 0,B = 1,X = 1 and k = I. 
Fig 10 : The phase diagram of the parameters depicting an attractor solution. The initial conditions 
chosen are v{0) = 0.2,^(0) = 0.3 (green); v{0) = 0.3,ti(0) = 0.3 (blue); v{0) = 0.4, u(0) = 0.3 (red); 
v{0) = 0.5,ii(0) = 0.3 (brown). Other parameters are fixed at a = 1,6 = 0.01, A = 0.3, = 1,A = 1 and 
K = 5. 




Fig. 11 Fig. 12 Fig. 13 Fig. 14 



Fig 11 : The deceleration parameter is plotted against the state parameter, e = 0.003 is considered. 
Fig 12 : The statefinder parameter r is plotted against the state parameter, e = 0.003 is considered. 
Fig 13 : The statefinder parameter s is plotted against the state parameter, e = 0.003 is considered. 
Fig 14 : The ratio of density parameters is shown against e-folding time. The initial conditions chosen are 
v{0) = 0.3, u{0) = 0.3. The other parameters are fixed at a = 1, ^ = 0, 6 = 0.5, B = 1, X = 1 and k = 1. 
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which gives the power law form of expansion of the universe. Like DGP brane here also in order to have an accelerated 
expansion of universe in RSII brane we must have < ^(_rs//) < |. Using this range of X{^iigii-^ in the equation 
qRSii _ _|_ we get the range of qi^^^^^ as — 1 < Qc^'^^^^ < 0. This is again consistent with an accelerated 

expansion of the universe as it was in DGP model. Like before here also we calculate the statefinder parameters {r, s} 
in order to get relevant information of DE and DM in the context of background geometry only without depending 
on the theory of gravity. The expressions of the statefinder pair (eq. in the RSII model can be obtained in the 

form 

riusn, = (l - (1 - iX.nsu,) • (53) 

and 

s{Rsii) = X^RSii) (54) 



4 Graphical Representation of the Phase plane Analysis 

Phase diagrams are drawn to determine the type of critical point in DGP and RSII models. We discuss the results 
obtained in the two models separately. 

DGP BRANE MODEL 

The dimensionless density parameters v and u are drawn in figures 1 and 2. From the figures we see that v 
decreases, and u increases during evolution of the universe. This shows that the density of DM decreases while the 
density of DE increases as the universe evolves. So this result is consistent with the well known idea of an energy 
dominated universe. The phase space diagram (figure 3) shows the attractor solution. The eigen values are calculated 
at the critical point and found to be (1.97625, — 1.09495). Hence the critical point is a saddle point. 

Figure 4 shows the variations of the deceleration parameter, q against Wmcg- From figure 4, it is evident that 
there is a gradual decrease in the deceleration parameter q, and finally in the late universe it attains negative values, 
which suggests that there should be an acceleration in the late universe. This result is in accordance with the various 
observational data of la supernovae and CMB data which suggests that the universe is undergoing an accelerated 
expansion of late. Figure 5 and 6 shows the variations of the statefinder parameters r and s against Wmcg- From these 
figures it can be seen that r tends towards 1 and s tends towards 0. Therefore it is evident that these results tends 
towards the KCDM model. Figure 7 shows the variation of the ratio of the density parameter v and u with time. 
From the figure it is evident that the ratio of the above parameters decreases with time. So it can be concluded that 
there is a relative decrease in matter density with respect to the energy density. This is again consistent with the 
notion of an energy dominated universe. 

RSII BRANE MODEL 

The dimensionless density parameters v and u in case of RSII brane model are drawn in figures 8 and 9. From the 
figures it is evident that at early stages of evolution of universe density of DM decreases steeply relative to the density 
of DE such that finally an energy dominated scenario is witnessed in the late universe. The phase space diagram (figure 
10) shows the attractor solution for the model. The eigen values are calculated at all the three critical points. At the 
first critical point the eigen values are found to be (402.877, 204.835). Hence the first critical point is an unstable node. 
At the second critical point the eigen values are found to be (1214.75, 565.097). Hence the second critical point is also 
an unstable node. At the third critical point the eigen values calculated are (—25.3905+2257.49*, —25.3905— 2257. 49i). 
Hence the third critical point is a stable focus. 



Figure 11 shows the behaviour of deceleration parameter q as it is plotted against the state parameter Wmcg- The 
figure shows that there is a gradual decrease in the deceleration parameter with time. Finally in the late universe the 
deceleration parameter q becomes negative, which confirms the fact that the universe is undergoing an accelerated 
expansion. In Figures 12 and 13, the statefinder parameters r and s are plotted against Wmcg- Finally in Figure 14 
the ratio of the density parameters v and u is plotted against time. From the figure it is found that just like DGP 
brane here also there is a relative decrease in matter density with respect to the energy density. This is consistent 
with the notion of an energy dominated universe. Henceforth we find a solution for the cosmic coincidence problem. 
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5 Study of Future Singularities 



It is a well known fact that any energy dominated model of the universe is destined to result in a future singularity. 
The study of dynamics of an accelerating universe in the presence of DE and DM is in fact incomplete without studying 
these singularities, which are the ultimate fate of the universe. MCG, being a very generalized form of DE can step 
into the phantom era with the EoS {w < —1) for definite values of the constants A and B. It is known that the universe 
dominated by phantom energy ends with a future singularity known as Big Rip (Caldwell, R. R., 2003), due to the 
violation of dominant energy condition (DEC). But this is not the only type of singularity possible. Nojiri et al (2005) 
studied the various types of singularities that can result from an phantom energy dominated universe. These possible 
singularities arc characterized by the growth of energy and curvature at the time of occurrence of the singularity. It 
is found that near the singularity quantum effects becomes very dominant which may alleviate or even prevent these 
singularities. So it is extremely necessary to study these singularities and classify them accordingly so that we can 
search for methods to eliminate them. The appearance of all four types of future singularities in coupled fluid dark 
energy, F{R) theory, modified Gauss-Bonnet gravity and modified F{R) Horava-Lifshitz gravity was demonstrated in 
(Nojiri, S. and Odintsov, S. D., 2011). The universal procedure for resolving such singularities that may lead to bad 
phenomenological consequences was proposed. Nojiri, S. and Odintsov, S. D., 2011 classified and studied four types 
of such singularities in detail. We analyze two possible future singularities in case of MCG in brane world scenarios. 
Type / ("Big Rip") : For t ^ tg , a— >-oo, p— >-oo and \p\ oo 
Type // ( "Sudden" ) : For t ^ ts , a ^ Us , p ^ ps and \p\ — >■ oo 

Here we have considered MCG as the DE in DGP and RS2 brane models. 1) Phantom era is realized for MCG 
when p < [^4^]°+''. From the MCG equation of state i.e p = Ap — ^ we see that when p — >• oo , |p| — >■ oo when 
a oo and t ^ tg- Therefore it is seen that Type I singularity is quite obvious for MCG. This is the well known Big 
Rip singularity. 

2) Suppose if p — >■ and ps ~ 0, then |p| —oo for t ^ tg and as Hence we get the type II singularity. This 
is called a sudden singularity. 

The stability of the critical points analyzed in the last section will be occurred much before than the singularities 
take place, so any direct impact of the singularities of the dynamical system will not be observed. The big rip and 
sudden singularity may be formed before the critical points to be formed. 

6 Concluding Remarks 

In this work, we have considered a combination of Modified Chaplygin gas in DGP and RSII brane gravity models. 
Our basic idea was to study the background dynamics of MCG in detail when it is incorporated in brane gravity. 
Dynamical system analysis had been carried out, critical points were found and the stability of the system around 
those critical points was tested for both DGP and RSII brane models. Graphical analysis was done to get an ex- 
plicit picture of the outcome of the work. In order to find a solution for the cosmic coincidence problem, a suitable 
interaction between DE and DM was considcircid. The dynamical system of equations characterizing the system was 
formed and a stable scaling solution was obtained. The numerical values of the density parameters for DE and DM 
almost coincided with each other (as seen in graphs 8 and 9). Hence this work can be considered to be a signifi- 
cant one as far as solution of cosmic coincidence problem is concerned. From the above analysis we conclude that the 
combination of MCG in brane gravity makes a perfect model for the expanding universe undergoing a late acceleration. 
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